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Abstract 

We consider cooperative reactions and we study the effects of the interaction strength 
among the system components on the reaction rate, hence realizing a connection between 
microscopic and macroscopic observables. Our approach is based on statistical mechanics 
models and it is developed analytically via mean-field techniques. First of all, we show 
that, when the coupling strength is set positive, the model is able to consistently recover all 
the various cooperative measures previously introduced, hence obtaining a single unifying 
framework. Furthermore, we introduce a criterion to discriminate between weak and strong 
cooperativity, based on a measure of "susceptibility". We also properly extend the model 
in order to account for multiple attachments phenomena: this is realized by incorporating 
within the model p-body interactions, whose non-trivial cooperative capability is investigated 
too. 

1 Introduction 

The phenomenon of cooperativity is widespread in biological and chemical sciences and it has 
been the focus of many theoretical and experimental investigations |16] , Cooperativity is a typ- 
ical "emergent" property, that directly links the description of a system at the single molecular, 
elementary level, with the macroscopic properties in complex macromolecules, cells and organ- 
isms. It is often a compelling task to exploit cooperative effects, such as amplification effects 
and high sensitivity to external parameters |4j|8j. 

A common feature of cooperative systems is the interaction among "active sites": This in- 
teraction can result, for example, in an increasing affinity for binding of substrate as more sites 
are occupied, referred to as positive cooperativity for binding; Conversely, if successive binding 
of substrate to active sites reduces the binding affinity and inhibits the occupation capability of 
another site, we speak of negative cooperativity for binding. 

Here we focus on the effects of cooperation on reaction rate curves: when measuring the 
progress of a reaction as a function of the concentration of substrate (e.g. the saturation of 
hemoglobin by oxygen) [T9] one finds that cooperation is typically associated to sigmoidal 



*Dipartimento di Fisica, Universita di Parma 
^Dipartimento di Fisica, Sapienza Universita di Roma 



1 



curves, in contrast with the hyperbolic Michaelis-Menten law |17j . holding when cooperativity 
is absent. Actually, the plethora of different phenomena where cooperativity effects have been 
observed in reaction rates is so spread that many different definitions have been suggested in 
the past, most of them being based on the measure of some kind of deviations from the classical 
Michaelis-Menten reaction kinetics |17) . While non-cooperativity is a well defined behavior in 
complex multi-sites binding, all the definitions of cooperativity still lack a unifying picture, 
despite sharing a common underlying mechanism. 

Here our aim is to suggest a general description of cooperative behavior for multi-sites systems 
by using a basic statistical mechanics approach, able to emphasize the role of interactions between 
sites and to recover all the different behaviors and definitions. In our model, binding sites are 
associated to two- values variables, indicating the binding/non binding state analogously to the 
classical Ising model [22j|6j. 

A first step in this direction has already been paved[22j by using linear spin chains for 
modeling nearest neighbors correlations among sites; here we move forward, by assuming a 
simple mean field interaction structure among active sites as this allows to understand complex 
behaviors of macromolecules. More precisely, minimizing a suitably introduced free energy, which 
depends on proper order parameters, we obtain an explicit expression of the reaction rate as a 
function of substrate concentration and of a possible coupling among binding sites. Within this 
framework, we can then clarify several general properties of cooperative binding, hence allowing 
a unified picture of cooperative behavior in reaction- rate curves. In particular, we show that the 
strength of the interaction among active sites discriminates in a simple way between Michaelis- 
Menten curves and cooperative behavior, the latter in agreement with several, well-known criteria 
(e.g. Hill, Koshland, range, etc.); furthermore we identify a critical interaction threshold where 
strong cooperativity sets in/out, also showing that a sigmoidal rate law is not necessary for 
cooperativity In our framework, also reaction rates curves with an almost discontinuous step 
in a defined range of concentration can be recovered. In particular, the behavior of the reaction 
curves around the discontinuous step can be given by an interesting representation in terms of 
"critical" parameters. 

Finally, by extending our model to multi-sites attachment (a reasonable assumption when 
binding occurs on macromolecules or on a homogeneous surface), we obtain reactions curves 
exhibiting different features (weak, strong, sigmoidal or discontinuous step), as observed experi- 
mentally and negative cooperativity may naturally arise. 

The paper is organized as follows: In section 2 we briefly summarize the principal definitions 
of cooperativity in chemical kinetics, and in section 3 we introduce our statistical mechanics 
framework. The latter is then throughout investigated as explained in section 4, where we also 
recover previous definitions of cooperativity. In section 5 we properly extend our model in order 
to account for multi-attachment too. Outlooks and conclusions then follow. 
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2 Measures of cooperativity 



As in elementary chemistry, cooperativity can be denned as the ability of ligand binding at one 
site on a macromolecule to influence ligand binding at a different site on the same macromolecule; 
according to whether the affinity of the binding sites for a ligand is increased or decreased upon 
the occupation of a neighboring site, we speak of positive or negative cooperativity respectively. 
A common example is provided by hemoglobin, which displays four binding sites whose affinity 
for oxygen is increased when the first oxygen molecule binds |22j. 

In order to study the existence and the extent of cooperativity, a convenient observable is 
the fractional saturation of binding sites at equilibrium 9 (9 = occupied sites/ total sites) for 
different values of free ligand concentration a0. Due to the number of possible forms and levels of 
cooperativity, a multitude of diagnostic tests have been introduced in literature (see for instance 
|17|). mainly based on the comparison to a non-cooperative behavior. The latter occurs when 
the binding sites (e.g. enzyme units) act independently and are identical in their activity; this 
behavior is observed for example in myoglobin in the presence of oxygen |22| . In such cases, the 
equilibrium reaction curve is well defined and exhibits the familiar hyperbolic Michaelis-Menten 
(MM) form 



where Km, called the MM constant, rules the affinity between components, and a accounts for 
the substrate concentration; note that 9m correctly saturates to 1 |16j . Hence, cooperativity is 
expected to be at work when the reaction rate does not follow MM kinetics. 

We now review different measures and, accordingly, definitions of cooperativity (see also 
|17ll5l[T8l[T]). neglecting those (e.g. Wong co-operativity |23|). which are restricted to particular 
forms for the reaction rates. We also omit those based on the shape of 9(a) near the origin, 
as this is not necessarily indicative of the behavior at intermediate substrate concentrations: 
Indeed, mixtures of positive and negative co-operativity, different degrees of co-operativity and 
the occurrence of tipping points can not be predicted by the behavior of 9(a) when a is close to 
zero. 

2.1 Koshland cooperativity 

The Koshland measure of cooperativity is one of the easiest to implement, being based on the 
coefficient k, defined as the ratio 



In chemical kinetics one is typically interested in finding the rate law governing a reaction, i.e. how the 
concentration of a given reactant or product varies in time. For elementary reactions (i.e. reactions with a single 
mechanistic step) the law of mass action is valid, and the reaction rate is proportional to reactant concentrations 
raised to a power defined by stoichiometric coefficients, but this is no longer true for complex reactions. However, 
if the substrate molecules react on binding sites to form a product, as is the case of enzymes, we expect the rate 
for product formation to be still proportional to the fraction of occupied sites 8. 




(1) 



3 



where ao.i and cto.g are denned as those values of ligand concentration such that 0(ao.i) = 0.1 
and, similarly, #(0:0.9) = 0.9: For the MM hyperbolic function 9m{&), one would get k = 81, 
independent of Km- Hence, one can notice that if k < 81 then 9(a) is expected to grow faster 
than 8m((*) and, according to Koshland, we speak of positive cooperativity, while if k > 81, then 
we have negative cooperativity. 

The index k has the advantage of being an absolute number, but the disadvantage of being 
sensitive to the scale of measurements and of depending only on ao.i and on ao.9) hence ignoring 
all information that can be derived from the shape of 9(a). 



2.2 Global dissociation quotient 

Another possible approach to measure the degree of cooperativity is in terms of deviations from 
the hyperbolic behavior by introducing the generalized equation |13| 

9(a) = ^— -. (3) 

v ; a + K(a) K ' 

In the absence of cooperativity K(a) is constant, viceversa one gets K(a) = a(l — 9)/9. Hence, 
the experimental value of K(a) can be interpreted as the ratio between free sites and occupied 
sites times a and it is often called "global dissociation quotient". The slope of this curve, sym- 
bolized by 7, i.e. 7 = dK(a)/da, is used as a quantitative measure of cooperativity: 7 equal, 
greater or smaller than zero would correspond to absence of cooperativity, negative and positive 
cooperativity respectively. Under this view, the phenomenon of cooperativity is equivalent to 
the change in the global dissociation quotient with ligand concentration; the larger the rate of 
change of K(a) and the larger the degree of cooperativity. 



2.3 Hill cooperativity 

An alternative way to quantify the deviation from an hyperbolic behavior is to use the so called 
Hill function, which provides another extension of the MM function and which reads 



a 



a h + K M ' 



(4) 



with Km and h constant. Indeed, Eq. (JJ|) is able to successfully fit most experimental data 
from cooperative systems, where the parameter h can be interpreted as the number of ligand 
molecules that bind to the macromolecule in a single step |15j . 

More generally, given experimental data for 9(a), one can define the Hill coefficient as 
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(5) 



which provides a measure of the extent of cooperativity: when n#(a) > 1 for all a > 0, we 
have positive cooperation; viceversa, if nn(a) < 1 for all a > we have negative cooperation. 
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For instance, for oxygen binding to hemoglobin (with N = 4 binding sites), the estimated 
Hill coefficient is about 2.8. Indeed, in biochemical processes rin is typically smaller than N, 
suggesting a mixed effect resulting from the interaction among sites 

It is worth noticing that the main difference between 7 and uh, is that the former represents 
the absolute change in global dissociation quotient per absolute change in ligand concentration, 
while the latter is related to the corresponding logarithmic changes. 

Finally, we underline that, given the experimental curve 6(a), one can straightforwardly get 
an estimate for njj, and this explains the popularity of such cooperativity measure, nonetheless 
nu constitutes a macroscopic property and it is not uniquely related to the difference in binding 
affinity due to the occupancy of next sites, which, conversely, is a microscopic property. 

2.4 Range cooperativity: weak or strong behavior 

The range measure is, again, based on a comparison between the reaction curve 0(a) and an 
"appropriate" hyperbolic MM curve, as summarized in the following (see also |17|): first, one 
specifies a value K > and builds the hyperbolic function Ok (a) = a /(a + K). Then, if 
6'(0) < 9' K (0) (0'(0) > 9' K (0)) and 6(a) and Ok (a) intersect at most once for all choices of 
K > 0, we say that 0(a) is strongly positive (negative) cooperative. If 0(a) and 0k(ol) intersect 
more than once, it is possible to formulate a concept of weak cooperativity (in the range sense): 
If 0'(O) < 0^(0) (0'(O) > 0^(0)) and 0(a) > K (a) (0(a) < K (a)) for large enough a, then 
0(a) is weakly positive (negative) cooperative. 

As proved in [I7J, for monotone increasing rate curves the range and Hill measures of coop- 
erativity are equivalent, more precisely, strong positive (negative) cooperativity with respect to 
the range measure is equivalent to a positive (negative) cooperativity according to Hill. However, 
differently from Hill measure, which requires the estimate of the slope for log 0(a) against log a, 
range measurements do not require interpolating the reaction curve from the data, and accord- 
ingly, the diagnosis of cooperativity is relatively simpler and also offers additional advantages 
|17| . In particular, the method is able to distinguish different degrees of cooperativity (weak and 
strong), even though such a distinction is still based on macroscopic evidences and cannot be 
directly ascribed to specific values or range of values in the microscopic parameters characterizing 
the phenomenon. 

3 Bridging chemical kinetics and statistical mechanics 

Let us consider a system (i.e. a big macromolecule, a homo-allosteric enzyme, or a catalyst 
surface) which has a set of iV interacting sites, numbered by an index i = 1, 2, . . . , N. Each site 
can bind one identical smaller molecule of a substrate, and we call a the concentration of free 
substrate molecules. When a site has a molecule bound on it, binding on all the other sites is 
enhanced (inhibited), which corresponds to a system with positive (negative) cooperativity. 

In analogy with the Ising model (see for instance |22|). we associate in complete generality to 
each binding site a dichotomic variable cij which takes the value +1 if the i-th site is occupied, 
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and —1 if it is unoccupied. A configuration of the molecule is then specified by the set of values 
01,02, •••,0'JV = W}- 

Once this bridge has been established, we can speculate on the way the binding sites interact 
with the substrate and/or among themselves: as the interactions are based on elettromagnetic 
exchange forces the third law of dynamics holds and a representation in terms of an Hamiltonian 
H]\[ is allowed^]. In this way we are able to ask for free energy equilibria of this Hamiltonian, 
namely to ask for the minimum energy and the maximum entropy principles: When the system is 
in equilibrium with a heat bath at a given temperature T, the probability of a configuration with 
energy Hn is proportional to the Boltzmann factor exp(—H^ /ksT), where kg is the Boltzmann 
constant and T is the temperature or, more generally, the level of "noise" experienced by the 
system. 

We first focus on the interaction between the substrate and the binding site which is expected 
to depend on both the substrate concentration and on the state of the binding site; we model 
such interaction by an "external field" h meant as a proper measure for the concentration of free 
ligand molecules. Then, for a non interacting model (whose potential cooperation features are 
obviously neglected) one can consider a microscopic interaction energy given by 

N 

H N ({a};h) = -hJ2^i- (6) 

i=l 

We notice that for h positive, Hn is negative and large in absolute value whenever there is a 
prevalence of occupied sites, while it is positive if most sites are unoccupied. On the contrary, 
when h is negative the smaller, negative, values for this one-body interaction correspond to a 
prevalence of unoccupied sites and vice versa. So we can think at h as the chemical potential 
for the binding of substrate molecules on sitea^|: When it is positive, molecules tend to bind 
to diminish energy, while when it is negative, bound molecules tend to leave occupied sites. 
If we assume that the free molecules are non interacting (which is strictly true when they are 
very diluted), the chemical potential can be expressed as the logarithm of the concentration of 
binding molecules and one can assume that the concentration is proportional to the ratio of the 
probabilities of having a site occupied with respect to that of having it empty |22| . In this simple 
case, being the sites non-interacting, the probability of each configuration is the product of the 
single independent probabilities of each site to be occupied and one finds 

a oc p(ai = +l)/p(<Ji = -1) = exp(+/i)/ exp(-/i) 

so that 

h = -\oga. (7) 

2 Strictly speaking the intrinsic symmetry of the exchanges implies symmetry in the couplings which ensures 
detailed balance and the latter ensures the equilibrium canonical framework to hold, ultimately an Hamiltonian 
representation of the phenomenon. 

3 We stress that our formulation of the problem allows to introduce concentrations as real variables (instead of 
a natural enumeration of the substrate elements), namely external fields, which has the mathematical advantage 
of retaining the canonical ensemble instead of the grand-canonical one. 
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Hence the limit h — > — oo corresponds to a vanishing concentration a, while when h — > +00 
the concentration a goes to infinity as well. It is straightforward to see that the reaction-rate 
associated to eq. ([6]) corresponds to the hyperbolic Michaelis-Menten curve [22]. 

Once we have successfully tested the statistical mechanics approach for the paradigmatic MM 
kinetics, we can proceed with our mapping and account for cooperativity just by extending the 
one-body theory previously outlined: To model the interaction between sites we simply introduce 
the full Curie- Weiss (CW) Hamiltonian [H] 



which, again, represents the microscopic energy of the system in a given configuration, so that 
the lower and the larger the probability of finding that configuration. Of course, J > 
means that there exists a two-body interaction between the binding sites, and the free energy is 
better minimized if sites "cooperate" by aligning consistently. 

The Hamiltonian (JHJ) was originally introduced to describe a system of atoms with magnetic 
moments: the spins <r's can point in two possible directions, h represents an homogeneous external 
magnetic field acting on each atom and spanning all real numbers, and J is a fixed number 
representing the two-bodies interaction strength (the exchange term), which for our model rules 
the cooperativity of the system. In fact, the first term in ([8]) is an interaction between the 
N(N — l)/2 couples of sites: if J is positive, this energy is smaller the larger the alignment 
displayed by the it's; In the following we focus on positive couplings Jij > 0. The first term 
ranges from its minimum value —J(N — l)/2, obtained in the extreme cases in which all sites 
are either all occupied or all empt}@, to its maximum value 0, obtained at half saturation of the 
sites. So this contribution to the energy is such that sites tend to behave coherently, and the 
constant parameter J fixes the scale of the microscopic binding energy. 

It is important to notice that our model corresponds to a so-called mean-field approximation, 
where each site interacts with all the remaining N — 1 with the same strength J: this assumption 
implies the minimal tractable oversimplification to investigate complex behaviors. Even more 
refined models describing for example the binding of molecules to semiflexible polymers, can be 
treated with a mean-field approximation in the limit of long-range interactions[12j. 

First of all, a stochastic approach like ours is meaningful only if the number of elements taken 
into account is large; in particular, our solution holds in the thermodynamic limit (i.e. N — > 00); 
this implies that the predictiveness of our model can be significant even for a finite size folded 
macromolecule, built up by, say, N ~ 1000 binding sites: The simplest way to realize a scaling 
able to retain the physical/chemical significance of the phenomenon under investigation^] is by 

4 Notice that this symmetry is actually broken in the presence of an external field, namely of a sufficiently large 
substrate concentration a which favors positive (occupied) states. 

5 In fact complex phenomena as a folding of a macromolecule or a sudden jump in the reaction rate curve could 
not be taken into account by modeling the system by a linear chain with only nearest neighbors interactions. The 
world of long range models -whose applicability we are testing here- seem to be more consistent, and eventually 
more realistic diluted systems may be studied in future works. 





i=l 
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assuming the number of neighbors for any site to scale with N and this is indeed consistent with 
the mean-field approach. In this way the statistical overlap among a global behavior and a local 
one is strong in the large N limit. Otherwise stated, the features of an infinite length molecule, 
whose interaction are properly rescaled within a mean field approach, can be still comparable 
with a realistic one. Furthermore, still among the benefits, by varying only the ratio J/T, our 
model provides a robust picture of the different regimes known for cooperative systems; it is also 
mathematically tractable and it accounts for a discontinuous rate law which linear models with 
nearest neighbors interactions cannot recover. 

From the other side, close to criticality, its prediction (which are unsensitive to the underlying 
topology) may barely represent the real phenomenon, quantitatively. However, as we will see, 
they can suggest an interesting measure to be performed on nearly discontinuous regimes, in the 
"critical" region. 

The model could be further extended by including, for instance, a given topology for binding 
sites [3l[2], in such a way that the effective number of neighbors per site is smaller than iV (and 
may even remain finite in the thermodynamic limit), or a coupling Jjj depending on the couple 
of sites considered, say, scaling with their distance; the resulting mathematics would be more 
complicated, without any real qualitative change. Thus, keeping the framework as simple as 
possible for the sake of clearness, we introduce the standard statistical mechanics definitions for 
the CW model. The normalized probability for a configuration {a} is then 

P N (W}; J, h, T) = Z N l exp {-H N ({a}; J, h)/k B T) (9) 

where the normalizing factor 

Z N (J, h,T) = ^ exp {-H({a}; J, h)/k B T) (10) 
M 

is termed the partition function. The number n{a} of occupied sites can be computed as 

N 1 

nW = ^-(l + ff! ), (11) 
i=i 

and the binding isotherm 9(a) is reconstructed from the average fraction of occupied sites (n)/N 
as a function of the concentration a 

6 = 1^ = ^11 n {°} p {<?) = \ + \m{J, h) (12) 
M 

where m(J,h) = ^AfTil/N represents the average magnetization per site in mean-field CW 
model, the brackets (.) account for the Boltzmann average (i.e. (.) = ^ CT . exp(— f3Hjy)/ exp(— /3iTjv)); 
the dependence on a is obtained simply by substituting h with (1/2) log a (see Eq. 7). Here m 
represents the average fraction of occupied (or unoccupied) sites with respect to half saturation. 
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In the limit of large N, the average value (n)/N can be obtained by minimizing, with respect to 
the parameter 9 G (0, 1), the effective free energy |14j 



F(9, a, T) = sup 



(--(28-l?--(l-9)\og(a)-Ts{9)) 



(13) 



The (only direct) independence of this variational prescription by the coupling constant (as the 
sup is taken just on 9) is the main strand we pave to obtain a unifying picture of the various 
catalogues of cooperativity in chemical kinetics: It holds whatever J. 

The first two terms at the r.h.s. of eq. (|13p stand for the internal energy, which corresponds to 
the Boltzmann average of the Hamiltonian H({a}, J, h) expressed as a function of 9 and a, while 



is the entropic term, whose weight is ruled by the temperature T. For small temperatures, 
the most likely configurations are those associated to small values of the internal energy. As 
we explained above, according to ligand concentration, Hjq is smaller either when there is a 
prevalence of occupied or empty sites, and it is maximum at half saturation. 

We also notice that the entropic term is connected to the logarithm of the number of config- 
urations associated to a given fraction 9 of occupied sites: It is maximized for fractions around 
1/2, which have a larger number of configurations associated, and minimized for 9 = 0, 9 = 1, 
corresponding to just one configuration and a vanishing entropy. The minimum of F is the opti- 
mal compromise between the minimization of the effective internal energy and the maximization 
of entropy. We expect that for large interactions J (or chemical potentials) the energy term in 
(|13p is the leading contribution, and the optimal fraction is ruled by this term, so that the sites 
tend to be in the same state and the binding isotherm displays a sigmoidal shape, while for small 
values of the interaction strength the leading term will be the entropic one, which at a fixed 
value of the chemical potential prefers disordered states, i.e. states obtained by a large number 
of configurations, and pushes the rate law towards a MM form: In this sense the temperature 
can be thought of as a noise because when it is high the system prefers to be in a disordered 
state. However, large or small temperatures are defined with respect to the interaction strength 
J and in the following we are going to fix the temperature equal to 1 and let the interaction 
strength vary to see all the possible regimes of binding isotherms arising from different values of 
J, as the latter rules cooperativity. 

The minimum condition for (I13p with respect to the order parameter 9 corresponds to the 
CW self-consistence equation [6J[14J 



and gives the average fraction of occupied sites corresponding to the equilibrium state for the 
system. As we will see later, this equation describes a second order phase transition when a = 1 
and J larger than the critical value J c = 1. Conversely, when J > J c and a = a c = 1, the 




(14) 




(15) 
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transition is first order and the average fraction 9 has a discontinuity. Equation (|15p holds 
rigorously just in the thermodynamical limit (N oo); one could also consider corrections of 
order (1/N) (see for instance |14j). but for large systems these are not important, the only relevant 
change is that at finite volume the discontinuous functions are mildly smoother (accordingly with 
the experimental counterparts |16j). 

Summarizing, the binding isotherm is a continuous function of a for every fixed J < 1, and it 
has a discontinuous jump in a = 1 for J > 1. The value J = 1 corresponds to a critical behavior 
where the derivative of 9 with respect to a becomes infinite at the reference concentration a = 1, 
and in the following we will see that this is connected to an infinite Hill coefficient!!. 



4 Obtaining chemical kinetics from statistical mechanics 

As discussed in the introduction, to determine the rate law one has to compute the average 
fraction of occupied sites as a function of the concentration of free ligand molecules a, and this 
dependence is encoded in the self-consistence equation (|15p . 

We consider separately the two cases < J < 1 and J > 1 for the interactions between sites, 
because, as stated, while in the former case 9 is everywhere continuous in a, in the latter it has 
a discontinuity in a = 1: it takes a value smaller than 1/2 when a — > 1~ and greater than 1/2 
when a — > I^Lj. In either cases, however, it is easy to check that 9(a) — > for a — > (which 
corresponds to the h —¥ — oo limit) and 9(a) 1 for a — » oo (h — > +oo). This means that the 
reaction rate vanishes when the substrate concentration vanishes and it saturates to one (all the 
binding sites being occupied) when the substrate concentration is large, as expected. 

When J — > 0, no cooperativity is expected (as the model reduces to a one-body theory) and, 
coherently, we recover the MM kinetics. 
In fact the rate (I15p can be equivalently expressed as 

aexp[2J(2g(q,J)-l)] 
J > ~ i + a eX p [2J(29(a, J) - 1)] [ ' 

which properly gives, for J = 

9(a,J)\ J=0 = -^. (17) 
1 + a 

Note that we do not lose generality when obtaining a/(l + a) instead of a/(K + a) at the r.h.s. 
of eq. (|17p : in fact we can rewrite the MM equation as 9(a) = K7,a/(1 + K^a) such that 
choosing Km ^ 1 is equivalent to shifting a — > a/ Km, which can be compensated by shifting 
even h (1/2) logo - (1/2) log 1^- 



6 An infinite Hill coefficient may seem unrealistic, however it is not an unavoidable feature of our modeling: 
in fact h scales with the connectivity of the underlying network of interactions and, while the latter diverges in 
this minimal fully connected representation, diluted mean fields can still work fine. However, as the mathematics 
involved becomes immediately very heavy we preferred to present the pure theory within the limitation of the 
high connectivity limit. 

7 In the magnetic counterpart this corresponds to the fact that below the critical temperature we have a 
ferromagnetic phase with a positive remanent magnetization when h — > + , and a negative one when h — > CP. 
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Moreover, from eq. (fT6|) . we get that when J > the rate for a given concentration is smaller 
than the corresponding one for a non interacting system when a < 1, and becomes greater when 
a > 1, in fact, the greater the interaction and the steeper the sigmoidal shape of the the rate. 
Equation [16] is plotted in Figure [1] versus a and for several values of J. 



4.1 Continuous rate-law 

As we said before, for J < 1, 9(a, J) is everywhere continuos^] and, when the concentration is 
equal to the reference concentration a = 1, one half of the sites is occupied in average. The 
shape of the kinetics law, however, depends on the interaction between sites. The derivativ^] 
of 9 with respect to a, which is strictly related to the Hill coefficient and, consequently, to the 
cooperativity of the system, can be computed from (fT5|) : 

89 _ 1 l-(2fl-l) 2 

8a 4a 1 - J [I - (26- l) 2 ] ' 1 ; 

It is always positive and finite for J < 1, meaning that 9 is an increasing function of a, as we 
expect. In the limit of low concentration we obtain 

= exp(-2J) (19) 

and the kinetics at very low concentration is governed by the two-bodies interaction J: the 
greater J and the flatter the rate law. Note that the strength of the cooperativity, J, appears 
in eq. (|19|) as an exponent, implicitly supporting the log-scale of the Hill coefficient in eq. (5). 
When J = 0, d a 9\ a= o = 1 and one properly recovers the same trend as that of the MM kinetics, 
which has a first order kinetics with the same coefficient for small concentrations. 



89 
da 



4.2 Weak and strong cooperativity 

The second derivative can be computed and expressed in terms of the first one: 



d 2 9 

da 2 



a da 



1 + 



29 - 1 



[1 — J[l — (29 — l) 2 ]) 2 



(20) 



For the magnetic model this is a paramagnetic phase, meaning that there is no remanent magnetization when 
the external field vanishes: m(a, J)| a =i = 0. 

9 Note that in the frame of the Curie- Weiss model this is strictly related to the generalized susceptibility 



X - 



dm(h) 
dh 



which measures the response of the system to a change in the field h. In fact, we have 

80 ldm(h(a)) 1 dh 1 
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Figure 1: Top: the figure shows the different phases for the system. For J < 1/4 the binding 
isotherms have no inflection point and the system is weakly (W) cooperative; when 1/4 < J < 1 
an inflection point arises, and the presence of a sigmoidal shape allows us to call this regime 
strongly (S) cooperative; both regions correspond to a continuous varying 0(a), while for J > 1 
6(a) becomes discontinuous, with a jump at the critical concentration a = l(blue line). Bottom: 
the different regimes for the binding isotherm obtained by varying the interaction strength J. 
For J = (blue line) the hyperbolic Michaelis-Menten law represents the isotherm for a non- 
interacting system; for J = 0.2 (green) the system is in a weakly cooperative regime; for J = 0.6 
(red) strong cooperativity manifests itself with the typical sigmoidal shape; J = 1 (light green) is 
the critical regime: the derivative in the inflection point which gives the Hill coefficient is infinite; 
J = 1.8 (purple) represents the discontinuous phase, with an extremely strong cooperativity. 
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When a ranges in (l,oo), m is positive and the second derivative is negative, so that 9 is a 
concave function of a in that range. For a = 1 we have d^ 2 f = — (l/4)/(l ~~ J), so that 9 is still 
concave there. For a € (0, 1) a numerical study of the second derivative clearly shows that it 
vanishes for a unique a* depending on J. This a* is near to a = 1 for values of the interaction 
close to J = 1, and it reduces with J towards a = 0. Interestingly, the inflection point does not 
disappear when J = 0, but for a greater value, namely J = 1/4. In fact, expanding m in the 
first order in a one finds 

|4 =-2(l-4J)ex P (-4J) (21) 

so, when 1/4 < J < 1, there is a unique inflection point a* depending on J, which separates 
the region where 9 is convex (small concentration), to the one where it is concave. For J = 1/4 
the inflection point corresponds to a* = and when J < 1/4 the rate function 9 is everywhere 
concave, tending, for J — >• 0, to the hyperbolic MM form (whose second derivative —2(1 + a)~ 3 
is always negative). Note that when J = the expression (|2ip gives, correctly, the MM value 
— 2. The absence of an inflection point in the region J € [0,1/4] allows us to define it as a 
weak cooperativity region, while we can call strong cooperativity the one arising in the interval 
J G [1/4,1]. 

These very simple definitions have the advantage of being directly related to the microscopic 
interaction, so that the experimental behavior of a system could allow one to reconstruct this 
interaction strength and interpret the rate law in terms of the very general mean-field model we 
introduced. 

Summarizing, we saw that the rate is continuous in a in the whole range < J < 1, and it 
has a unique inflection point in the region 1/4 < J < 1, and no inflection points for J < 1/4 
(see Figure [T]). As a sigmoidal rate has necessarily an inflection point, we talk about strong 
cooperativity in the former case and weak cooperativity in the latter. 



4.3 Strong cooperativity and discontinuous rate-law 

The different phases for the binding isotherm are shown in figure Figure [TJ When J < 1 we can 
consider two regions, a weak cooperative one, where the rate law is hyperbolic, and a cooperative 
one with an inflection point for the rate law that grows gradually with the interaction strength. 
When J > 1 (corresponding in the original Ising model to the "ferromagnetic" phase) the rate 
law is still increasing with a, and the expressions (|19H2ip remain valid for a/1. In this point 
the rate function is discontinuous and the jump is given by 0+(J) — 9-{J), where 

0±{J) = lim 9(a,J). 

a— ►l* 

These two limits depend on J: they are both equal to 1/2 for J = 1, when the curve is still 
continuous, and their difference increases smoothly with the square root of J — 1 when J > 1. 
This means that, starting from vanishing concentration, the system has less sites occupied, for a 
given q, than the corresponding non interacting one, until the concentration reaches the reference 
value. Here, it is sufficient to increase infinitesimally the number of free molecules to obtain a 
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large filling (depending on J). After that value, the number of occupied sites is always greater 
than the corresponding value for MM. Note that, in principle, if the concentration varies slowly 
one could observe met ast ability, with a curve which continues growing continuously up to values 
of a > 1. The entire out of equilibrium features of the model are ruled out in this treatment 
as we deal with equilibrium statistical mechanics, however -as a second step- the bridge could 
be extended in that direction. If J — > oo this discontinuity increases, while its derivative in 
zero vanishes, so that in the large volume limit we obtain a step function. This corresponds 
to a chemical kinetics where no binding site is occupied until the concentration has reached 
the critical value a = 1, and when this value has been reached all sites are occupied. This 
kind of discontinuous behavior can be observed, for example, in the binding isotherms of small 
surfactants onto a polymer gel |20| . 

When J — > 1 a second order phase transition appears. This indicates that the correlation 
between binding sites becomes stronger and the typical trend of thermodynamical observables is 
a power law |14j . If the statistical mechanics picture we suggest applies to the reaction rates, then 
this power law behavior becomes very interesting, as it should be related, as in the statistical 
mechanics counterpart, to very general features the systems, due to the concept of universality 
|14| . Let us consider in details the behavior for J — > 1. As we said before, the discontinuity for 
J > 1 is given by #+(</) — #_( J), whose dependence from J near the critical point (a = 1, J = 1) 
can be expressed as 

+ (J)-0_(J)«(J-l)V2 

while on the critical isotherm (i.e. for J = 1) around the critical concentration a = 1, mean-field 
theory predicts 

0(a,l)-^(a-l) V3 . 

In this regime one can also predict the behavior of the a dependence of the Hill coefficient defined 
in ([5]), which, when a — 1, diverges as 

n H (a) « (a- 1) _2/3 . 

Moreover, we know that when J — > 1^ the susceptibility, and so the derivative of 9 respect to a, 
diverges as 

x\ a =i ~ ii - j\~ • 

As hinted in Sec. 3, the exponents of these power laws are only valid in the limiting case of 
a very large number of interacting neighbors per site, and we expect that the right exponents 
should depend on the real dimension of the space in which they are embedded, as in the cor- 
responding statistical mechanics models. However these scalings, in particular those related to 
the reaction rate around the discontinuity as a function of a, suggest a new interpretation of 
almost discontinuous reaction curves and could represent an interesting measurements to test 
our predictions. Indeed, statistical mechanics suggests that power-law behaviors in the vicinity 
of the discontinuous step are signatures of a critical phenomena: If in particular the detected 
exponents do not differ that much from the ones predicted by the CW model, then this suggests 
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that the mean field picture, despite its simplicity, is able to capture the essence of the cooperative 
phenomenon. 



4.4 Cooperativity through the Hill coefficient 

As explained in Sec. 12. 3| the usual way to define in a quantitative manner the cooperativity of 
a system is by the Hill coefficient, i.e. the slope at the symmetric point 9 = 1/2 



1 
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9(l-9)dloga 
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6=1/2 



log a 



6=1/2 



(23) 



If binding on different sites is an independent process, one simply finds uh = 1, while in the 
extremum case in which sites are either all empty or all occupied nu = N. We call a system 
cooperative (non cooperative) if uh > 1 (tih = 1)> while the cooperativity is said to be negative, 
meaning that binding is reduced if there are occupied sites, for njj < 1. So this gives a lower 
bound for the number of interacting sites, and it is possible to see that it is related to the variance 
of the mean number of occupied sites. 

The Hill coefficient for our general model depends, as expected, on the interaction J; in 
particular for J < 1 we have 



log [0/(1-0)] 
log a 



4 — 

Oct 



1/(1 - J) 



(24) 



a=l 



Being the derivative of for a = 1, the Hill coefficient is finite (and greater than one) for 
J < 1 and it diverges for J — > 1~ when the discontinuity appears. We observe that the last 
equation relates nu to the typical scale of interaction energy. This is a macroscopic measure of 
cooperativity which is directly associated to the microscopic interactions among sites. 



4.5 Cooperativity through the global dissociation quotient measure 

Among the useful tools to describe cooperativity, we recall the global dissociation quotient K(a), 
whose derivative 7 is expected to be different from zero when some form of cooperation occurs; in 
particular, 7 should be negative when there is a positive cooperation. Figure [2] shows some plots 
of 7(a) for different values of J: consistently with the definition of K(a), (see !2.2p . in the region 
of a < 1, the stronger the interaction, the smaller 7, properly indicating the deviation from the 
inconstant MM curve. As the derivative of with respect to a appears in the definition of 7, in 



10 When there is not such a symmetry, one should more properly consider the quantity 

39 



tih = 4 



d log a 



(22) 



max, mm 



that is the slope of #(loga) calculated at the inflection point, where it has an extremum. In fact for an unsym- 
metrical system the slope at 9 = 1/2 is not generally an extremum slope and should not be taken as a suitable 
cooperativity index. 
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correspondence of the critical value J = 1, 7(a) diverges for a = 1. In fact this is the (critical) 
region where the correlations between sites, and so the cooperativity, are expected to diverge. 

For sake of clarity in the figure, we did not show a plot of 7 for J > 1, however in this case 
the curve is below the others shown for a < 1 and it is not defined for a = 1: from this point 
on, it is practically vanishing. 




Figure 2: Top. The figure shows different regimes for the derivative of the global dissociation 
quotient 7(a) = dK/da obtained by varying the interaction strength J. For J = 0.2 (green) 
7 remains quite close to zero, indicating that K is nearly constant in a, as one expects for 
a non-cooperative system; for J = 0.6 (red) the global dissociation quotient has a stronger 
dependence on a: we know that in this case the binding isotherm has a sigmoidal shape and 
strong cooperativity appears; when J = l(cyan) we are in the critical regime, where the derivative 
of 9, and so 7, diverges at a = 1. Bottom. The Koshland coefficient k = ao.9/ a o.i obtained by 
our theory as a function of the two-bodies interaction J. It is a decreasing hyperbolic function, 
defined for values of J below J = 1.5, which tends to be lower when cooperatively is larger and 
the binding isotherm takes a steeper sigmoidal form. 

4.6 Cooperativity through the Koshland measure 

Lastly, we want to recover even the measure of cooperativity introduced by Koshland (see sec. 
2.1): To satisfy this task, we plotted the previously defined Koshland coefficient as a function of 
the two-bodies interaction coupling J (see Figure [2]) , as the latter is the only relevant tunable 
parameter to explore cooperativity at work. 
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The coefficient is defined only for J < 1.5; in fact for J > 1 9(a) is a discontinuous function and 
when J > 1.5 the equations 9(a) = 0.1 and 9(a) = 0.9 cease to have a solution. 
As shown in Figure [2] (bottom), we find that it is a decreasing function of J, which takes the 
value k = 81 for non-interacting systems (J = 0), and k ~ 1 when J = 1.4: Coherently with 
his scenario, we note that 9(a = 1) = 1/2 for every J < 1, when J > the binding isotherm is 
below the MM isotherm for a < 1 (such that ao.i is larger than in MM theory) and above MM 
curve for a > 1 (such that cto.9 is smaller than in MM theory). 

5 Multiple interacting systems 

When considering the binding of systems like long chain molecules, possibly on a homogeneous 
surface |10) , it may be possible that the interaction involves more than two hosting sites per time 
and this can in principle affect the global "cooperative capability" of the system. Thus, in the 
following, we want to extend our model by accounting for multiple interactions, involving more 
than two elements together, in particular three and four bodies interactions. 

The simplest p-body interaction among N sites which can be either occupied or empty (cr$ = 
+1 or —1 respectively) can be expressed by the Hamiltonian (12511 

oo . N 

H N ({a}, J,h) = -^2 2NP- 1 Jp 2 a h ...a iN -h^ai. (25) 

p=2 l<i 1 <...<i p <N i=l 

The interaction strength J is assumed to be positive (for J = we recover the standard MM 
kinetics again) and the combinatorial factor before the summation makes the energy extensive 
and accounts for a sensible p — » oo limit |7j. 

It is easy to see that the energy (the mean value of the Hamiltonian) scales as 

p-l 

(H N ) (x N(m$(m)), $(m) = ^ c k m k + h, (26) 

k=i 

where with we meant the coefficient of the Taylor series implicitly defined by eq. (|26|) . and 
m = 29-l. 

The argument in the self-consistent equation m = tanh(#) is then built up (in complete 
generality) by all the terms of the Taylor series (whose convergence is tackle by the p\ term in 
the numerator). 

However, as \m\ < 1 we expect that only the first terms do matters in the thermal average 
and the main contribution is given by the p = 2, which has been studied in details so far. 

However, let us highlight here another advantage in performing the statistical mechanics 
approach: As the Hamiltonian 1251 represents a sum of Hamiltonians (each for each p), which, 
in principle, can have coupling strengths very different (i.e. the distribution of J p can range 
over several orders of magnitude), in this case solving the global problem as a whole can be 
prohibitive, while here we can consider each Hamiltonian contribution by itself and focus only 
on one term every time (namely p = 2, 3, 4, ....): 
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This procedure, which we are going to perform, allow to isolate the different contribution 
(and weight in the global behavior) of each independent binding capability represented by each 
independent p-term in the sum. 

Note that the a's acting on a single site could tend to keep it in the state —1 or +1 in 
accordance to their product. So, for example, for p = 4 one could have three sites, acting on 
the i-ih site, in a configuration, say (+1, +1, -1), which favors the state <jj = —1, even if in 
the whole there are more sites in the state a = +1. The same local field is obtained when the 
configuration is (-1,-1,-1), i.e. if all the three sites are empty. 

The energetic behavior of this extension is then deeply different from the previous case: 
Usually one deals with linear forces, which are generated by quadratic potentials (i.e. p = 2) 
such that the Maxwell-Boltzmann probability distribution is Gaussian accordingly with what 
naively expected by a simple Central Limit Theorem argument and critical behavior (at h = 0) 
arises to confirm this picture [6]. 

However if a more complex scenario appears then a violation to this linear framework is 
expected: for instance for p = 3 it is straightforward to check that even a positive J may have 
subtle anti-cooperation features: in fact it is immediate to check that the energy would prefer the 
orientation of the spins +1, +1, +1 but also +1,-1,-1 (this can be thought of as a competitive 
feature of the multi-attachment that naturally introduce negative cooperativity in the process 
under investigation). 

Interestingly, this extension still shares with the simplest p = 2 case the same entropy: In 
fact, as in the two-body case, for this long range interacting system the energy can be easily 
expressed as a function of the parameter 9 describing the fraction of occupied sites, and of the 
concentration a 

N- l H p {6,a) = - J -{29 -l) p -\ log(a)(20 - 1). (27) 

while the entropy per site is exactly the same as in the p = 2 case. Given the effective free energy 
(still keeping T = 1 for the sake of simplicity) 

F(6, J, a) = sup ( - ?-{28 -lf-\ log(a)(20 - 1) - s(9)) , (28) 
the minimum condition with respect to the order parameter 9 reads off as 

0(J, a )~\ = \ tanh - l) p ~ 1 + \ Ma)) . (29) 

This corresponds to the equilibrium state for the system, and the average fraction of occupied 
sites will be given by the solution of this equation. Once again, this equation, as in the p = 2 
case, is strictly valid when the number of sites is large. 

Again, we expect that for large interactions J (or chemical potentials) the energy term in 
(|28p is the leading one, and the sites tends to be in the same state (this corresponds to a large 
magnetization) , while for small values of the interaction strength the leading term is the entropic 
one, which prefers disordered states, i.e. states where the sites do not see each other. 
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5.1 Case p = 3 




Figure 3: p = 3. Top: the figure shows the phases of a p = 3 system. For J3 < 0.50 the system 
is weakly (W) cooperative, with no inflections; for 0.50 < J3 < 0.86 = J% c it is strongly (S) 
cooperative, having at least one inflection point; the blue line shows the ^-dependent critical 
concentration at which, for J3 > J^ c , a discontinuity in 6(a) appears, passing from a relatively 
small 0(1) to a larger one (II) as the concentration increases (at fixed J3). Bottom: binding 
isotherms 6(a) for different interaction strengths J3. For J3 = 0.4 (blue line) the system is 
weakly cooperative, with a hyperbolic form; for J3 = 0.8 (green) it has two inflections, the first 
for a = 1/2 and the second for a larger a; for J3 = 1.1 (red) the isotherm has a strong inflection 
for a < 1 and a discontinuity for a value of the concentration between 1 and 1.5, corresponding 
to a transition to a more ordered state; for larger J3, as the cyan curve for J3 = 1.5 shows, the 
discontinuity is shifted towards a < 1. Note that all the curves pass through the point (1, 1/2). 

In this case the energy (|27j) is an odd function of m = 26 — 1 at fixed a 

N- l H(6, a) = -J 3 (26 - l) 3 /2 - i log(a)(20 " ( 30 ) 
Its global minimum, corresponding to a fraction of occupied sites 6 and satisfying 

= \ + \ tanh (I J 3 (26 - l) 2 + ~ log(a)) , (31) 

is located at an average number of occupied sites smaller than 1/2, when the concentration is 
very small. On the contrary, when the concentration is around a = 1 (so that the chemical 
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potential for binding is small) and J3 2> 1 the global minimum of energy corresponds to a large 
number of occupied sites. There is an intermediate region where H (9, a) can have two local 
minima for 9 smaller and larger than 1/2 respectively, but the relative strength of the chemical 
potential log (a) versus the interaction J3 tells us which is the global one, corresponding to the 
equilibrium state. It is easy to see that the two minima (located respectively around 9 = 1/4: 
and 9 = 1) correspond to the same energy on the line a = exp(— 3J3/4) (this is the critical 
line where there is a two phases coexistence in the limit of large J3 or log (a) and of negligible 
entropy contribution (T — > 0)). However the meaningful function to minimize is the whole free 
energy and in general one has to take into account the entropy term, which tends to favor an 
equal number of occupied or empty sites and to disfavor deviation from this. 

So we can identify three regions in the (J^,a) plane: When J3 is very small, the leading 
interaction is guided by the chemical potential log(a), so that the sign of 9 — 1/2 is the same 
as the former and it vanishes for a = 1. This corresponds to a weakly (W) cooperative system. 
Growing J3, the system begins to feel mutual interactions, and the average number of occupied 
sites as a function of a has an inflection point for a* > 1, increasing quickly before that point 
and more slowly after (see for instance the case J3 = 0.8 in Figure This is due to the interplay 
between entropy, which pulls the system towards 9 = 1/2, and the energy, which prefers 9 = 
or 9 = 1 when the chemical potential is sufficiently large (Figure EJ. 

As in the previous p = 2 case, we can perform a numerical study of the second derivative 
with respect to a of the binding isotherm, in order to discriminate between weak and strong 
cooperativity, to understand when an inflection point, and consequently a sigmoidal shape, ap- 
pears: While for p = 2 we found that this happens for J = 1/4 in a = so that this exact value 
for J can be inferred by developing the second derivative for small values of a, here this is no 
longer true, as we find that the first inflection point is obtained for J3 = 0.5, in correspondence 
of a = 1.5. This is related to the fact that the minimum of the free energy passes quite abruptly 
(however continuously and it is ultimately due to a violation of the quadratic shape for the ener- 
getic term) from values of 9 not much larger than one half to larger values. We can then identify 
J3 = 0.5 as the value separating a weak cooperativity regime from a strong cooperativity one. 
When J3 is larger than this value, another zero in the second derivative appears, so that we have 
two inflection points. Above the critical value J3 = 0.9, instead, we have a unique inflection 
point below a = 1. 

For a critical interaction strength of J3 = 0.86, the system has a first order transition so 
that 9 is continuous in a below a critical concentration a*(J3) > and it grows abruptly above 
this threshold. Increasing J3, this critical concentration becomes lower up to negative values 
at certain J3. In this case, due to the strong interaction, a concentration smaller than a = 1 
in sufficient to have this discontinuous transition between a negative value of 9 — 1/2 and a 
positive one. When J3, and so the cooperativity is large, the critical line tends to be located at 
a = — exp(3J3/4). 
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5.2 Case p = 4 



As for p=2, in this case the energy is symmetric with respect to 9 = 1/2 when the concentration 
is equal to one, so that also the free energy features this symmetry 

F(0, a) = sup ( - 2 J 4 (29 - l) 4 - - log(a)(20 - 1) - s(9)) (32) 
\e\ v 2 ' 

and thermodynamic stability requires its minimization (w.r.t. 9) as 

6 = \ + \ tanh(2 ,h(29 - l) 3 + ^ log(a)). (33) 

When the interaction among sites is small, again, we have a binding fraction ruled by the con- 
centration of free ligands and vanishing as the latter goes to zero. For J4 = J c \ = 0.69 however, 
the system has a first order phase transition so that, coming from small concentrations, there 
is a value a c (J4) < 1 for which 9 changes abruptly to a larger value, which is smaller than 
1/2. Then 9(a) varies continuously in the interval a c (J 4 ) < a < a c (J±)\ passing through the 
point (1, 1/2), and there is a new transition with a discontinuity at an inverse value of a c (Ji)-, 
Oi c {Ji)' = a c (J^) . In fact for this range of interaction strengths J4, and for concentrations 
around a c (J 4 ), the free energy has two local minima coming from the well matched competition 
between the energetic and the entropic terms, the former preferring small, "ordered" 9, while 
the latter has always its maximum at 9 = 1/2. The concentration tells us which is the global 
minimum and the two minima are equal at a c . The same happens around the symmetrical value 
a c (J 4 )'. 

The critical concentration a c tends to one with the interaction strength growing, and for 
J4 = J C 2 = 1-37 the system becomes ferromagnetic. In this case the mutual interaction is the 
overwhelming force driving the system, and, as in the previous p = 2 case, the binding isotherm 
has a discontinuity when passing from a < 1 to a > 1, changing abruptly from a large negative 
value (depending on J4) to a positive, opposite value (see Figured]). It is interesting to note 
that the values of a c and a' c depend on J4 and even for p = 3 we find a critical concentration 
depending on J4: because this dependence is absent in the p = 2 case, an experimental J4- 
dependent critical concentration may indicate that a multiple interaction effect is acting (this 
can be seen for example in |20j). 

As for p = 2,3, we can determine exactly the value J4 = 1/24 as the one for which a first 
inflection point appears, so that the binding isotherm passes from a weak to a strong cooperativity 
region. In fact, as for the p = 2 case (see eq. l2Tj) by gaining the second derivative of the binding 
isotherm for small a one can see that the leading term is proportional to (1 — 24J4). The isotherm 
has a unique inflection point until J4 = 0.45, then a second point appears at a = 2.4, which for 
larger J4 splits in two points, so that for 0.45 < J4 < J c \ we have three inflections. The ones 
corresponding to the smaller and to the larger values of a disappear when J c \ is reached and, 
instead of these inflections, the isotherm has two discontinuities. 

A discontinuous behavior which could be explained in term of two and multisite interactions 
between hosting sites has been observed in the binding isotherm of long chain alkyl sulphates 
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and sulphonates to the protein bovin serum albumin, and in the adsorption isotherm of alkylam- 
monium chlorides chains on the biotite surface |9j . Typically, if there are multi-site interactions 
(with p > 2) it is very likely that also two-bodies interactions are present, and one should con- 
sider the several possible interactions with different strengths J p in the energy term, while for 
the sake of clearness we considered them separately. 

Before concluding some remarks are in order: 
As we have shown, the phenomenology pertaining to systems with p-body interactions depends 
sensitively on p and this allows to infer information about the properties of the system under 
study, starting from its reaction rate. 

On the other hand, we also notice that when looking at the most likely configuration, that 
is the most likely value of 9 for a given parameter set (p,J,a), we find that it exists and it 
is unique: if we, for the sake of simplicity, focus just on even p > 2 values, we can see that 
the amount of solutions is constant, namely, independent by the amount and complexity of the 
cooperating/anti-cooperating binding sites. More precisely, the self-consistent relation in Eq. 28 
may allow for one or two distinct solutions, but only one is a global minimum for the free energy: 
the minima of the free energy do not scale, but simply shift, with p, while the global minimum 
is always unique. This argument may be applied to the problem concerning the folding of (long) 
proteins, whose secondary and tertiary structure is essentially always the same, despite a large 
number of multi-attachment is in principle possible |21j . 

6 Conclusions and perspectives 

In this paper we pursued the strand early paved by J. J. Thompson |22j in modeling chemical 
kinetics with statistical mechanics techniques: highlighting the emergent properties in this field, 
namely cooperation among binding sites, we relaxed the original assumption of dealing with 
geometric models (i.e. linear chains) in favor of mean field models as the latter may properly 
account also other experimental findings (as discontinuous jumps or folding) and benefit of a 
broad tunable scenarios able to recover essentially all the definition of cooperativity introduced 
in literature. 

After a streamlined introduction to the modern definitions of "cooperativity" in chemical 
kinetics, we developed a one-to-one mapping between chemical variables (binding sites, sub- 
strate, etc) and the statistical mechanics counterparts. Then, we applied statistical mechanics 
prescriptions, that is the variational principle (minimization of the free energy) and we found an 
expression for the reaction rate 9 (representing the number of occupied sites) as a function of the 
substrate concentration a. In doing so the interaction constant J between binding sites is treated 
as a free (microscopic) parameter. We show that by tuning J we are able to recover a plethora 
of different behaviors, ranging from non-cooperative (J = 0), weak cooperativity (J < 1/4) and 
strong cooperativity ( J > 1/4). Our results are also shown to be consistent with the number of 
definitions of cooperativity (e.g. Hill, Koshland, or range), previously introduced. Furthermore, 
we found a universal criterion to distinguish between weak and strong cooperativity, that is, by 
looking at the inflections of the second derivative of the free energy. Indeed, we found that the 
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existence of cooperativity (i.e. is a positive interaction constant) is not sufficient for the existence 
of a flex in the reaction rate. Otherwise stated, a non-sigmoidal law does not necessarily mean 
that no cooperativitiy is at work. Such situations are here referred to as cooperativity of weak 
degree. 

There are several directions along which this work can be extended: for example including a 
well defined topology for binding sites, which would account to more realistic interaction scalings 
among the binding sites. Also, non-strictly-positive couplings would realize a spin-glass system 
which, as well known from statistical mechanics, significantly enriches the phenomenology, to be 
possibly mapped in terms of reaction kinetics. 

Acknowledgements 

The research belongs to the strategy of exploration funded by the FIRB project RBFR08EKEV 
which is acknowledged. 

References 

[1] Luis Acerenza and Eduardo Mizraji. Cooperativity: a unified view. Biochimica et Biophysica 
Acta, 1339(1):155-166, 1997. 

[2] Elena Agliari and Adriano Barra. A hebbian approach to complex-network generation. EPL, 
94(1):10002, 2011. 

[3] Elena Agliari, Adriano Barra, and Federico Camboni. Criticality in diluted ferromagnets. 
J. Stat. Mech., page P10003, 2008. 

[4] Daniel J. Amit. Modeling brain function. Cambridge University Press, 1987. 

[5] W. Bardslesy. Shapes of velocity curves in multiunit enzyme kinetic systems. Biochem. J., 
153(101):97-115, 1976. 

[6] Adriano Barra. The mean field ising model trough interpolating techniques. J. Stat. Phys., 
132(5) :787-809, 2008. 

[7] Adriano Barra. Notes on ferromagnetic p-spin and rem. Math. Meth. Appl. Sci., 32(7):783- 
797, 2009. 

[8] William Bialek and Sima Setayeshfar. Proc. Natl. Acad. Sci., 102(29): 10040-10045, 2005. 

[9] J. M. Cases and Frederic Villieras. Thermodynamic model of ionic and non-ionic surfactant 
adsorption-abstraction on heterogeneous surfaces. Langmuir, 8(5)4251-1254, 1992. 

[10] Avik P. Chatterjee and Pinaki Gupta-Bhaya. Multisite interactions in the lattice gas model: 
shapes of binding isotherm of macromolecules. Molecular Physics, 89(4): 1173-1 179, 1996. 



23 



[11] P. F. Cook and W. W. Cleland. Kinetics of Enzyme Mechanics. Oxford Press, 2007. 

[12] Haim Diamant and David Andelman. Binding of molecules to dna and other semifiexible 
polymers. Phys. Rev. E, 61:6740-6749, 2000. 

[13] J. T. Edsall and J. Wyman. Biophysical Chemistry. Academic Press, New York, 1958. 

[14] R. S. Ellis. Entropy, large deviations and Statistical Mechanics. Springer Press, 2005. 

[15] Archibald V. Hill. The possible effects of the aggregation of the molecules of haemoglobin 
on its oxygen dissociation curve. J. PhisioL, 40:4-7, 1910. 

[16] Terrell L. Hill. Cooperativity Theory in Biochemistry. Springer, Berlin, 1985. 

[17] Sam Karlin. Measures of enzyme co-operativity. J. Theor. Biol, 87(l):33-54, 1980. 

[18] Sam Karlin and Ron Kenett. Shapes of velocity curves in multiunit enzyme kinetic systems. 
J. Theor. Biol, 52(l-2):97-115, 1980. 

[19] A. Marangoni. Enzyme Kinetics: A Modern Approach. Wiley Press, 2002. 

[20] Y. Murase, T. Onda, K. Tsujii, and T. Tanaka. Discontinuous binding of surfactants to a 
polymer gel resulting from a volume phase transition. Macromolecules, 32:8589-8594, 1999. 

[21] A. Pagnani, G. Parisi, and F. Ricci-Tersenghi. Glassy transition in a disordered model for 
the rna secondary structure. Phys. Rev. Lett, 84(9):2026-2029, 2000. 

[22] Colin J. Thompson. Mathematical Statistical Mechanics. Princeton University Press, 1972. 

[23] J. Wong. Kinetics of Enzyme Mechanics. Academic Press, New York, London, 1975. 



24 



Figure 4: Top: phases for the p = 4 model. Binding is continuous for J4 < J c % = 0.69, 
and it passes from a weak (W) to a strong (S) - both continuous - regime at J4 = 1/24; for 
J c i < J4 < J C 2 = 1-37 isotherms have two discontinuities, one for a critical a c < 1 and the 
other one for a a' c = l/a c > 1 (corresponding to a critical line, so that we can identify three 
(I, II, III) continuous regions at a fixed J4; when J4 > J C 2 however, there is only one large 
discontinuity at a = 1, and the behavior is similar to the one observed in the J4 > 1 phase 
for p = 2 (see Figure [T]). Bottom: binding isotherms 9(a) of a 4-bodies interacting system, for 
different values of the interaction J4. The binding isotherm for J4 = 0.04 (the blue line) has no 
inflections and represents a practically non cooperative system; when J4 = 0.5 (green) the curve 
has an inflection point for a small concentration; for J4 > J c \ two inflection points appear and 
the isotherms develop two discontinuities corresponding respectively to the concentrations oi\ 
and 02 such that a\a% = 1, see the curves for J4 = 0.8(red) and J4 = l.l(cyan); for interactions 
larger than J C 2, the system has just one discontinuity for a = 1 and the binding isotherm tends 
to a step function, as for the purple line corresponding to J4 = 1.5. 
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